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This is a status report of the evaluation of the three-loop corrections to the static QCD potential of a heavy 
quark and an antiquark. The families of Feynman integrals that appear in the evaluation are described. To 
reduce any integral of the families to master integrals we solve integration-by-parts relations by the algorithm 
called FIRE. To evaluate the corresponding master integrals we apply the Mellin-Barnes technique. First results 
are presented: the coefficients of nf and n^, where n; is the number of light quarks. 
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1. Introduction 

The QCD potential between a static quark and 
its antiquark can be cast in the form 
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where the renormalization scale of a^ is set to q . 
The one-loop contribution ai is known since al- 
most 30 years |1I2] and also the two-loop term has 
already been computed end of the nineties |3|4|5j . 
Furthermore logarithmic contributions are known 
at three- and four- loop level [6|7| . Explicit results 
are nicely summarized in the recent review [5]. 
The non- logarithmic third-order term, as, is still 
unknown. 

2. Reduction to master integrals 

Any Feynman integral that contributes to a^ 
can be mapped to one of the three graphs shown 
in Fig. [1] where solid lines stand for usual mass- 
less propagators of the form l/(— p^ — iO)"' and 
wavy lines stand for linear propagators l/{~v ■ 
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Figure 1. Families of Feynman integrals needed 
for the calculation of 03. The solid lines corre- 
spond to relativistic massless propagators and the 
zigzag lines represents static propagators. 



p — SiiO)"'% with V ■ p — poi Si = ±1 and integer 
indices at . In case the latter type of propagators 
is absent the integrals reduce to usual massless 
two-point functions which can be treated with the 
help of MINCER [9]. Note, however, that the pres- 
ence of the static lines significantly increases the 
complexity of the problem. 

In general the integrals involve up to fifteen 
propagators (including an irreducible numerator). 
In order to simplify the reduction problem we ap- 



ply in a first step partial fraction identities to 
arrive at various subfamilies of integrals with at 
most three linear propagators. Thus any result- 
ing integral is labeled by twelve indices one of 
which corresponds to an irreducible numerator 
and three indices correspond to linear propaga- 
tors. 

Let us in the following describe the Feynman 
integrals that are generated by the n; contribu- 
tion. Altogether we have to consider about 70 000 
integrals (allowing for a general QCD gauge pa- 
rameter ^). As far as the "ladder" and "non- 
planar" diagrams are concerned we have to deal 
with the type of integrals shown in Fig. [2] where 
the linear propagators appear in the following 
form: If the loop momenta, k, I and r, in Fig. [2] 
are chosen as the momenta of the three upper 
lines, then the first diagram appears in two ways: 
either with the product 

(-U • fc - iOy^'i-v ■ I - ioy^^i-v ■ r - iO)-"", 

or with the product 

{-vk + iOy<'{-v ■ I - iOy^^i-v ■ r - ioy^K 

The second diagram appears with similar prop- 
agators where the momenta {fc,/,r} (in the first 
variant with — iO in three places) are replaced by 
{k,k — ?,r}. The third diagram in Fig. [2] cor- 
responds to {fc, Z, Z — r} and the fourth one to 

{k,l,T}. 

The integrals which correspond to the "Mer- 
cedes" graph are shown in Fig. [31 If we choose the 
loop momenta fc, I and r as the momenta of the 
three lower lines these five diagrams appear with 
linear propagators of the form {—v ■ k — iO) with 












Figure 2. Ladder and non-planar diagrams con- 
tributing to the ni part of 03 . 



Figure 3. Mercedes diagrams contributing to the 
ni part of a^. 



momenta {fc, k~r, I}, {r,k — l,r — I}, {fc, r,k — I}, 
{fc, r, I}, {fc, r — I, I}, respectively. 

The next step is a reduction of all the nine types 
of these Feynman integrals to master integrals by 
solving integration- by-parts relations (lOj. To do 
this we apply the algorithm called FIRE (Feynman 
Integral REduction) |llll2ll3ll4irM B] which is 
based on an extension of the classical Buchberger 
algorithm to construct Grobner bases (see, e.g., 
Ref. [II]). 

Similarly to other approaches, we work in a 
given sector, i.e. a domain of integer indices at 
where some indices are positive and the rest of the 
indices are non-positive. The aim is to express 
any integral from the sector in terms of master 
integrals of this sector and integrals from lower 
sectors, where at least one more index is non- 
positive. It turns out that in the higher sectors 
(with a small number of non-positive indices) the 
corresponding s-basis |12I13I14I15J (a kind of a 
Grobner basis) can be constructed easily (and, in 
most cases, even automatically). 

In the opposite situation where a lot of non- 
positive indices occur, s-bases are constructed not 



so easily. Usually there is the possibility to ex- 
plicitly perform an integration over some loop 
momentum for general value of e with results 
in terms of gamma functions. A straightforward 
way to do this leads to multiple summations and 
turns out to be impractical. An advanced strat- 
egy within FIRE is to use s-bases for some re- 
gions of indices corresponding to a subintegral 
over such a loop momentum in order to reduce 
these indices to their boundary values. Then it 
is sufficient to use explicit integration formulae 
only for the boundary values. Integrals which are 
obtained from initial integrals by an explicit inte- 
gration over a loop momentum in terms of gamma 
function usually involve a propagator with an an- 
alytic regularization by the shift e or 2e. After 
this integration we obtain a two-loop reduction 
problem with seven indices which is then solved 
by FIRE. 

Finally, after using Grobner bases in higher sec- 
tors and an explicit integration in lower sectors, 
it is still necessary to solve the reduction prob- 
lem in a relatively small number of intermediate 
sectors. In these cases we turn to Laporta's algo- 
rithm [18I19J implemented as part of FIRE. 

To reduce by FIRE all the integrals contributing 
to the ni part of 03 it took around ten days on a 
2.3 GHZ Opteron computer with 8 GB operative 
memory. The most complicated integral in this 
reduction is F{\, ...,1,1, —4, 1, 0) corresponding 
the nonplanar graph in Fig. [5] where the index of 
the irreducible numerator is zero and one of the 
static propagators is raised to the power —4. All 
other indices are equal to 1. 

3. Evaluating master integrals 

After using FIRE we know the master integrals. 
The number of the master integrals appearing in 
the ni part of 03 is around one hundred. For 
their evaluation we used the Mellin-Barnes (MB) 
technique which is based on a replacement of a 
sum of terms raised to some power by their prod- 
ucts in some powers, at the cost of introducing 
additional integrations over contours of a com- 
plex plane. Then one takes explicitly all other 
integrations (over loop momenta and/or over al- 
pha/Feynman parameters) and is left with a mul- 



tiple MB integral. 

For planar diagrams, experience shows that a 
minimal number of MB integrations is achieved if 
one introduces them loop by loop, i.e. one derives 
a MB representation for a one-loop subintegral, 
inserts it into a higher two- loop integral, etc. 

Consider, for example, the dimensionally reg- 
ularized Feynman integral corresponding to the 
first graph in Fig. [2] with the linear propagator 
{—vk+iO)'^°-^ . Let us denote it by F{ai, . . . , an). 
A straightforward implementation of the loop- by- 
loop strategy leads to a six-fold MB representa- 
tion which reads 
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r(-6 -I- ai,...,8 + gg, 10, 11/2 + 3g + ^2,3) 
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r(6 — ai,2,3,4,5,7,8 — ag, 10, 11/2 — 3e — ^2,3,4) 

r(a6 - 2:4)r(a5 - Z5) 
r(a2, 3, 4,5,7,8 + Qio,ii/2 + 2e-A + ^1,2,3,4) 
r(6 - 02,3,4,5,7,8,10 - aii/2 -3e - Z2 + zq) 
r(2 - Qe - a9/2 - s + zi,4)r(a5 -I- 23,4 - Z5) 
r(a2, 3,4,8 + aii/2 + e -2 + 22,5,6) 
xr(a4 + Z2,5)r(aio/2 + zi - Z6)r(aii/2 -|- zq) 
xT{2 - 04,8 - aii/2 - £ - Z5,6) 
xr(2 - 03,4 - aii/2 - £ - Z2 + ze) 
xr(2 - 05,7 - aio/2 - £ - zi,4 -t- Z5,6) 
xr(a3,4,8 + aii/2 + £-2 + Z2,5,6), (2) 

where 03,4,8,11 = 03 -I- 04 -f a8 + an, zi,2,3 = 
2i + Z2 -l- Z3, etc. By definition, any integra- 
tion contour over Zi should go to the right (left) 
of poles of gamma functions with +z-dependence 
( — z-dependence) . 
There are two strategies for resolving the singu- 



larities in e in MB integrals suggested in Refs. PUI 
m] (see also Chapter 4 of ;22|23j ). The second one 
was formulated algorithmically [24I25J . and the 
corresponding public code MB . m ^5" has become 
by now a standard way to evaluate MB integrals 
in an expansion in e. Using MB.m and evaluating 
the resulting finite MB integrals by corollaries of 
Barnes lemmas we have obtained, for example, 
the following result for one our master integrals: 
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In contrast to a similar Feynman integral with the 
linear propagator l/(— w • k — iO) (see Ref. ^5] ) 
one needs, at a first step of the resolution of the 
singularities in e, an auxiliary analytic rcgulariza- 
tion. In fact, the result in Eq. ^ and the corre- 
sponding one for the case of l/(— w • k — iO)°^ (see 
Ref. [55]) differ by a term which is nothing but 
a residue taken when performing the initial ana- 
lytical continuation to the auxiliary parameter of 
analytic regularization. 

The above integral is an example where prod- 
ucts of the type l/(x -I- iO)°l/(a; — iQ)^ appear, 
with X = —V ■ fc, —V ■ /, etc. If a and b are pos- 
itive integers, these products are ill-defined. For 
example, they do appear in diagrams resulting 
from iterations of the Coulomb potential. How- 
ever, for integrals in Eq. ^, and many other in- 
tegrals that we encounter, in these products at 
least one of the exponents is regularized by an 
amount proportional to e and/or z (i.e. a vari- 
able of the MB integration). As a consequence 
they are well-defined, in the sense of analytic con- 
tinuation. More explicitly, we have 
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sin(7r(a -I- b)) 

for general complex numbers a and b. 

All the master integrals that we encounter are 
expressed in terms of C(*)i 1—2,3,4,5,6, powers of 
In 2, Lij{l/2),j ~ 4,5,6, and the constant sg 



(see, e.g., PTl) up to transcendentality level six. 
With the assumption about the presence of these 
numbers in the results, one can apply the PSLQ 
algorithm '28] in situations where explicit integra- 
tion over MB parameters by corollaries of Barnes 
lemmas is no longer possible. If this is just a one- 
fold MB integral then it is possible to obtain a 
numerical result with more than 300 digits which 
are sufficient for the PSLQ algorithm. For twofold 
MB integrals, one can hope to obtain an accuracy 
of 50 digits depending on the complexity of the 
integrand. This accuracy can also be sufficient for 
a successful application of the PSLQ algorithm, 
at least when a given integral possesses homoge- 
neous transcendentality. Unfortunately, sufficient 
criteria about homogeneous transcendentality are 
not known at the moment so that this property 
can be seen only experimentally, by considering 
lower terms of the expansion in e. 

4. Preliminary results and perspectives 

The three-loop correction to the static quark 
potential can conveniently be parameterized in 
the form 



as 



(3) 3 , (2) 2 , (1) , (0) 

flg Ui + a^ Hi + a^ rii + a^ , 
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where n; denotes the number of massless quarks. 
Using the techniques described above we evalu- 
ated the coefficients a^ and 03 which read 
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The nf contribution together with the CaTp part 
has already been presented in Ref. [26] . 

The coefficient Og ' is also reachable. Only four 
constants contributing to the higher-order expan- 
sion in £ of some master integrals are not known 
analytically at the moment. 
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